Chiral Extrapolation, Renormalization, and the Viability of the Quark Model 
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The relationship of the quark model to the known chiral properties of QCD is a longstanding 
problem in the interpretation of low energy QCD. In particular, how can the pion be viewed as 
both a collective Goldstone boson quasiparticle and as a valence quark antiquark bound state where 
universal hyperfine interactions govern spin splittings in the same way as in the heavy quark systems. 
We address this issue in a simplified model which; however, reproduces all features of QCD relevant 
' to this problem. A comparison of the many-body solution to our model and the constituent quark 

model demonstrates that the quark model is sufficiently flexible to describe meson hyperfine splitting 
provided proper renormalization conditions and correct degrees of freedom are employed consistently. 
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One of the fundamental symmetries of the strong interaction is chiral symmetry (xS). This symmetry arises because 
the bare masses of the light u and d quarks are small with respect to hadronic scales (or, m u> d « ^-qcd)- The 
nonzero quark masses imply that chiral symmetry is explicitly broken. Chiral symmetry is also broken spontaneously 
by the QCD vacuum. This fact is represented by the nonzero value of the quark condensate which plays the role of 
. an order parameter. 

The phenomena associated with chiral symmetry breaking, for example the spectral properties of quarks in the 
dressed vacuum or the formation of Goldstone modes, are among the most important and most extensively studied 
aspects of strong QCD. As a consequence of chiral symmetry, the axial current is (partially) conserved and the 
interactions of low momentum pions with hadrons are weak. This forms the basis of chiral perturbation theory and, 
more generally, enables one to formulate an effective approach to low energy QCD with dynamics being dominated 
Qh| by weakly interacting Goldstone bosons and nucleons [jj. At high temperatures or densities the chiral invariance of 
the QCD vacuum is expected to be restored. The thermodynamics of this phenomenon plays an important role in 
determining signatures of the quark-gluon plasma in relativistic heavy ion collisions || . Finally, xS* is also relevant 
for lattice gauge simulations. In particular, chiral symmetry predicts a highly nonlinear dependence of a number of 
observables on the light quark mass. Since current lattice QCD simulations are limited to light quarks with masses of 
. i-h ! the order of 200 MeV the interpretation of lattice results requires input from chiral symmetry-based phenomenology. 

It is clear that chiral symmetry plays a central role in determining the structure and interactions of low energy 
hadrons ||. 

The constituent quark model (CQM) has historically developed in parallel with chiral theory and has been exten- 
sively used as a simple alternative for the study of hadronic phenomenology. Even though it is loosely connected 
to QCD, the quark model has had a large number of remarkable successes Q. The CQM is based on the idea that 
strong interactions lead to massive quasiparticles (constituent quarks) and that hadronic structure is dominated by 
the interactions between valence constituent quarks. One argues that constituent quarks are the effective degrees of 
freedom arising after dynamical \S breaking due to bare quark interactions with the quark sea in the chiral nonin- 
variant vacuum In the CQM, however, properties of the constituent quarks, e.g. masses and magnetic moments, 
are treated as free parameters making the approach effectively insensitive to the underlying chiral structure of QCD. 
This implies, in particular, that CQM pions lose their nature as Goldstone bosons and are not much different from, 
say, p mesons. The mass splitting between the pion (a spin-0 constituent QQ state) and the p meson (a spin-1 state) is 
attributed to a residual hyperfine interaction. The hyperfine interaction is typically associated with the nonrelativistic 
reduction of single gluon exchange between constituent quarks and yields a vector-pseudoscalar splitting proportional 
to a s /(m q m q ). 

Since the notion of gluon exchange is intrinsically perturbative, it is likely a useful concept in heavy quark systems 
where the large quark mass guarantees its applicability. However its relevance to the light quark sector is unclear. 
Indeed, one may expect that the ir-p mass splitting is largely driven by the underlying x$ rather than by the 
hyperfine interaction. It is therefore surprising that when the splitting between the flavored quark-antiquark vector 
and pseudoscalar mesons is plotted against the mass of one of the constituent quarks (the other constituent quark 
mass is held fixed), it does indeed behave like l/m con typical to a hyperfine interaction all the way down to the light 
quark sector ||. This is shown by the squares in Fig. 1. The definition of the CQM masses depends on details of the 
model, here we use typical values m u — raid — 330 MeV, m s = 550 MeV, m c = 1600 MeV, and mt = 4980 MeV. 
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FIG. 1. Flavored pseudoscalar and vector meson mass splittings as a function of the renormalized heavy quark mass (solid, 
mq = ulr) and constituent quark mass (dashed, tuq = m CO n) together with lattice results [12] (triangles, mg = ran) and 
experiment [13] (squares, mq = m con ). The dot-dash curve represents the hyperfine splitting obtained when the hyperfine 
interaction is neglected. The renormalized light quark mass is fixed at 5 MeV. 

The question we address here is whether the apparent universality of the effective interactions of the constituent 
quark model is consistent with the dictates of chiral symmetry in the light meson sector. To address this issue it is 
natural to work within the framework of Hamiltonian Coulomb gauge QCD. This enables us to map the degrees of 
freedom and the global symmetries of QCD to those of the constituent quark model and thus effectively "derive" the 
quark model from QCD. The Coulomb gauge QCD Hamiltonian can be formally written as Jt],|| H = H can (A)+6H(A). 
Here the first term is the canonical Hamiltonian which contains the Coulomb interaction Hcouiomb regularized at a 
scale A. The cutoff A is necessary since products of field operators at the same point lead to divergences. Modifying 
the structure of the canonical Hamiltonian by the cutoff requires counterterms, SH which are to be adjusted so that 
physical observables are A-independent. We note that in the chiral limit A is the only scale in the Hamiltonian and 
that light meson masses will therefore be proportional to the cutoff. 

Even though one can in principle choose an arbitrary value for the cutoff, there is typically an optimal choice, A = An 
which depends on the basis and approximation scheme used to diagonalize the Hamiltonian. Loosely speaking, as 
A increases beyond An, the strength of the various interaction terms decrease but at the same time the available 
phase space for partons to interact increases. On the other hand, for small cutoffs the phase space decreases but the 
remaining Fock sectors interact strongly, making it difficult to determine the effective degrees of freedom. Since our 
goal is to understand the quark model from QCD, we choose A^ to match to the quark model. The quark model 
scale is the scale where all QCD dynamics are given in terms of effective interactions between valence constituent 
quarks (that is if the quark model is to be a rigorous consequence of QCD). Thus, in particular, since soft gluons are 
not present in the Fock space at the quark model scale there are no contributions to hyperfine interactions from soft 
gluon intermediate states. 

We now examine the appropriate An for heavy- heavy and heavy- light systems. Since the heavy- heavy case is 
analogous to Coulombic bound states, it is natural to choose An so that it scales with tur = m(A — An) - the 
heavy quark mass renormalized at the scale Ar. This constrains the average momenta of heavy quarks to lie below 
tur. Furthermore, since interactions between heavy quarks and transverse gluons vanish in the nonrelativistic limit, 
gluons and heavy quarks decouple and the dominant interaction between heavy quarks is due to the nonabelian 
instantaneous Coulomb potential. As discussed in detail in Ref. ||, in nonrelativistic QCD (NRQCD) the proper 
scaling of the quark-gluon interaction is obtained if Ar = ctRrriR where an = cx(Ar) = (p)/mR. Here (p) denotes the 
typical momentum of a heavy quark in the bound state. The scaling (p) ~ aRiriR and E — 2m,R ~ a R m,R for average 
momenta and energies is typical for a Coulombic bound state. 

We now consider the heavy-light case. In this case the dynamics of the "brown muck" of gluons and light quarks is 
still subject to the full complexity of nonperturbative strong interactions. A heavy-light quark system with a single 
heavy quark is no longer Coulombic and thus the average separation between the quarks increases causing them to 
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probe the long range confinement potential. Quark model phenomenology flCfl dictates that this potential be linear. 
Thus the relevant quantities are the scales of the linear potential b and the kinetic energy, which is characterized 
by the reduced mass [i. Dimensional analysis then indicates that the average momenta in the bound state scale as 
(p) = (6/i) 1 / 3 . Since the reduced mass is very nearly equal to the light constituent quark mass, and this mass and \/b 
are of the order of Aqcd, we simply use A R = (p) Aqcd for the heavy-light case. 

Proper determination of the renormalization scale is the first step toward "derivation" of the quark model from 
QCD. The second, is to implement dynamical chiral symmetry breaking which will generate masses for the light bare 
quarks thereby creating the constituent quarks of the CQM. The constituent Fock space then provides an efficient 
basis for diagonalizing the Hamiltonian. 

In this letter we do not aim at a rigorous derivation of the meson spectrum. Our goal is to illustrate how the two 
steps discussed above enable one to derive the constituent quark picture while maintaining correct chiral properties. 
We will therefore approximate H and use a simple contact interaction instead of the full nonabelian Coulomb kernel. 
This is sufficient as long as the interaction causes the typical momenta in the bound state to scale as (p) = a R m R for 
heavy-heavy and (p) = Aqcd for heavy- light mesons. 

Hcouiomb -» H C (A) = J dx [^ t (x)T^(x)V t (x)T^(x)] A (1) 

Here the subscript A indicates that the operators are to be point-split or smeared over a distance ~ 1/A | and c c (A) 
is the coupling which replaces a (A). The counterterms SH contain relevant, marginal, and irrelevant operators. In 
the chiral limit, chiral symmetry prevents the occurrence of a relevant operator in the quark sector. For finite quark 
masses, the relevant operator is absorbed into the definition of the bare quark mass. The effect of transverse gluons 
eliminated by the cutoff show up through contact operators of dimension six or greater. There are a number of such 
terms; to illustrate the effect of hyperfine interactions we consider the dominant spin-dependent term, 

H h (A) = ^ J dx [^(x)Ta^(x)V + (x)TaV(x)] A , (2) 

where Ch(A) is proportional a(A). Thus the full model Hamiltonian is given by H = J i/j'[—ia ■ V + f3m(A)]tfj + 
H c {A)+H h (A). 

As stated above, generating spontaneous chiral symmetry breaking is crucial to examining the interplay of the 
hyperfine interaction and the Goldstone boson nature of the pion. As a result, we employ the BCS ansatz to construct 
an approximation to the chirally noninvariant vacuum. For the above Hamiltonian this leads to the following mass 
gap equation which determines the effective constituent quark mass 

^con / ~ / a \ - /a\\ / 2 



m (A) + -^(c c (A)-c fe (A)) / g'dq (3) 

yi J v m con + q 



where c c = Cfc c /2tt 2 , Ch = 3CFCh/2n 2 , and the single quasiparticle energies are given by -E'(q) = \/ m con + Q 2 - These 
resemble energies of the constituent quarks if m con ~ 300 MeV for the light quarks. 

We wish to carefully distinguish the three types of quark masses which have appeared. They are the bare cutoff- 
dependent quark mass, to (A) which is defined as the parameter appearing in the mass term of the canonical Hamil- 
tonian, the renormalized quark mass, m R , given by m(A R ), and the constituent quark mass, m coni which is the 
nonperturbative mass obtained by solving the gap equation. The cutoff dependence of the gap equation and m con are 
discussed later. 

The RPA or Bethe-Salpeter equation for the meson bound state, M, is then given by |ll[] 

(M\[H, Q M ]\BCS) = (E M - E BCS )(M\Qi\BCS), (4) 

where Q^ M is defined in a standard way in terms of the positive and negative energy wave functions, Q M = 
^apBliDp — ip'pDpBa with B and D being the quasiparticle operators. In the simple approximation to 

H used here, Eq. (Q) is an algebraic equation for the bound state masses Em = Em(A, m(A), c c (A), c^(A)). Requiring 
the Em to be A-independcnt is used to determine the cutoff dependence of the couplings. 

In the following we will concentrate on the pseudoscalar (ps) and vector (v) open flavor mesons. As discussed above, 
in the case of unequal quark masses, the renormalization constant should be fixed at A^ « Aqcd- Here we use the 
numerical value 420 MeV because this corresponds to \fb. The average, renormalized light quark mass, mji = m(A^) 
is set to 5 MeV, and the two renormalized couplings cr^Ji = c c ,h(An) are determined by fitting the it and p meson 
masses, E„ tP = E^^ p (A Rl m Rl c R:C ,c Rth ). 
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Solving the gap equation yields a light constituent quark mass of 380 MeV. For a heavy-light meson, the solid line 
in Fig. 1 shows the dependence of the hyperfine mass splitting A/^y s = E v — E ps as a function of the renormalized 
heavy quark mass tuq = vcir with the renormalized light quark mass fixed at 5 MeV. As ttiq increases the splitting 
falls off as l/m<g, however, as the heavy quark mass approaches the light quark limit the slope changes, reflecting the 

emergence of chiral symmetry: E ps oc ttiq 2 . Our predictions reproduce the lattice results (triangles in Fig. 1) |l2j ] 
available for large quark masses. The lattice results are given for f s ■ a as a function of ttir ■ a where a is the lattice 
spacing. Lattice calculations match our results if one takes a" 1 = 1.45 GeV which is very close to that of Ref. Jl4) . 
It is nontrivial that our predictions match the lattice over a large range of quark masses since the coupling constants 
were fixed in the chiral limit. We have fit the ansatz Mhf s — a + b\/ mn , + d + cniR to the data of Fig. 1. This 
ansatz has been employed in Ref. |l5| to fit their lattice data and was found to be accurate for low quark masses. We 
similarly find it to be a very accurate representation of the mass splitting up to a quark mass of roughly 1 GeV. 

The dashed-dotted line corresponds to the mass splitting when the hyperfine interaction is set to zero and the 
strength of the confining term is increased so that the light constituent quark mass remains unchanged. It follows 
that even in the chiral limit roughly 80% of the n — p mass splitting is due to the presence of the hyperfine interaction. 
This ratio is 92% and 99% for K* — K and D* — D respectively. We stress that the numerical value of this ratio 
depends on the details of the confining interaction. It is possible, for example, to obtain the correct ir — p mass splitting 
without a hyperfine interaction by adjusting the strength of the confining potential. The large ttlq behavior of the 
vector-pseudoscalar mass difference would however not be properly described, ie., it would decrease with ttiq more 
rapidly than 1/toq. In this case it also found that the confining interaction alone leads to predictions for the quark 
condensate and pion decay constant which are too small. For example, typical predictions for the condensate are 
roughly -(95 MeV) 3 |r|; which should be compared with the phenomenological value of -(250 MeV) 3 . In this model 
we obtain (qq)(An) = —(200 MeV) 3 . The improvement is due to the hyperfine term, which has not been accounted 
for in previous calculations. 

We now examine whether it is possible for the quark model to mimic the chiral behavior of the hyperfine splitting 
shown in Fig. 1. Our hypothesis is that the freedom to define quark masses in the quark model is sufficient to 
reproduce the splitting. We address this by plotting the hyperfine splitting as a function of the constituent quark 
mass derived from Eq. 3 (this is shown as a dashed line in Fig. 1). The curve reproduces the observed splitting for 
p — 7T, K* — K, D* — D and B* — B, shown as squares in the figure. We conclude that it is possible for a CQM 
to mimic the effects of chiral symmetry breaking. This is true because chiral symmetry breaking creates massive 
quasiparticles which may be used as effective degrees of freedom in model building. Furthermore, the hyperfine 1/toq 
behavior which is valid for heavy quarks continues to be valid for lighter constituent quarks. 

We emphasize that the solid line in Fig. 1 is a prediction which may be compared with lattice data (triangles). 
The only adjustable parameters, the two renormalized potential strengths, are fixed by the pion and rho masses when 
using 5 MeV for the renormalized light quark mass. 

We now consider hyperfine splittings in quarkonium systems as a function of the quark mass. In order to effectively 
mimic a Coulombic bound state for large mR and a confined state for small the renormalization scale is set to 
A_r = CRTOft + (m^b) 1 / 3 + 6 1 / 2 . In contrast with the heavy- light open flavor system, the hyperfine splitting does 
not vanish with increasing heavy quark mass ttlr but becomes proportional to c R mn, as expected for a Coulombic 
system. We also considered the case where the renormalization scale corresponds to purely "confined" systems, i.e., 
Ar = (m^fc) 1 / 3 + 6 1 / 2 . This leads to Mhf s oc A R /mn oc m R ^, a behavior which is different from that expected in the 
constituent quark model (where Mhf s oc A R /m 2 R oc m R ). The difference can be traced to the presence of negative 
energy solutions in the RPA equation. Elimination of the negative energy solutions leads to an effective hyperfine 
interaction which is proportional to c R h A 2 R /mR, while alone leads to CR.hA R /m R . Thus it is clear that correctly 
choosing the renormalization scale is crucial to establishing contact with quark model phenomenology. 

We conclude with a discussion of the renormalization procedure. In the simple model studied here there are three 
renormalization scale-dependent quantities: m(A), c c (A), and c^(A). We have used the pseudoscalar and vector 
meson masses to determine the (nonperturbative) renormalization scale dependence of these parameters. The third 
renormalization condition is obtained by fixing m con (or alternatively m(Afi) as done above for the light quarks) 
and by taking m con to be A-independent. For large A/Ar Eq. (^) implies that the bare quark mass behaves as 
to(A) oc m 3 on ln(A/m co „)/A 2 . Furthermore, at large cutoff the couplings saturate to constants. Thus asymptotic 
freedom is not recovered in our model. This is because the contact central interaction leads to matrix elements of order 
c c (A). In contrast, the full nonabelian Coulomb potential yields matrix elements of order c c (A)ln(A). Asymptotic 
freedom is obtained because renormalization then requires that c c (A) oc 1/ ln(A) for large A. Instead of choosing m con 
to be a physical parameter one may use the condition d[m(A)(qq)(A)]/dA — 0. This leads to a running constituent 
mass m con = m con (A) which vanishes for large A. 

The use of the contact interaction rather than the full Coulomb interaction leads to several other differences in 
phenomenology. For example, the constituent quark mass becomes a function of the quasiparticle momentum when 
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the Coulomb potential is employed. Single particle energies are enhanced (even divergent) due to the confining nature 
of the nonabelian Coulomb potential. Thus, in a realistic calculation one cannot use the constituent quark energy 
(or constituent mass) to constrain parameters of the interaction. Instead one may, for example, use the condition on 
the quark condensate given above. Finally, the contact interaction leads to a single bound state per channel which 
certainly precludes analysis of the corresponding quark model excited states. While all of these points affect details 
of our computation, they do not change the main conclusion. 

Spontaneous chiral symmetry breaking causes a pseudoscalar-vector meson mass splitting, but the hypcrfinc inter- 
action is required to obtain the necessary I/thq behavior for both heavy and light constituent quark masses. The 
interaction between the constituent quarks need not retain a memory of the underlying chiral symmetry and is de- 
scribed well by an effective hyperfine interaction. Thus Fig. 1 represents a direct validation of the main assumption 
of the naive quark model. 

We have shown how quark model phenomenology may be derived from a simple model of QCD. Elimination of high 
momentum components from quark-transverse-gluon coupling leads to short range hyperfine interactions. By choosing 
a renormalization scale which matches the quark model we neglect contributions from long range interactions (due to 
the exchange of low energy gluons) to the hyperfine interaction - as is consistent with the quark model. Then, studying 
the quark mass dependence of the pseudoscalar-vector mass splitting we are able to show that interactions between 
constituent quarks derived from QCD indeed follow that of the naive quark model while respecting chiral symmetry. As 
a result we have shown how the heavy quark mass limit extrapolates to the chiral limit and have illustrated the interplay 
between hyperfine interactions and chiral dynamics. The potential of utilizing standard many-body techniques in 
applications to QCD in the Coulomb gauge has been discussed by the authors and many others (|Jl^Jl^]; however, 
neither the role of renormalization in constructing the Hamiltonian and building the constituent quark representation 
nor a QCD-based demonstration of the applicability of the effective quark model interactions to both heavy and light 
quarks has been addressed in this context before. Even though the analysis presented here has used a simplified central 
potential, the key results are independent of this choice. This is because the simplified model captures all features 
(i.e., the presence of hyperfine interactions, correct momentum scales and dynamical chiral symmetry breaking) of 
QCD which are relevant. A study employing the full nonabelian potential will be presented elsewhere. 
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was supported in part by DOE under contracts DE-FG02-00ER41135, DE-AC05-84ER40150 (ES), and DE-FG02- 
87ER40365 (AS). 
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